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ABSTRACT: The constraints imposed by symmetry on the open channels 
of dipositronium has been studied, and the symmetry-adapted lowest open 
channel of each quantum state has been identified. Based on this study, the 
existence of two more + bound states has been theoretically confirmed, and 
a + narrow resonance has been predicted. A variational calculation has 
been performed to evaluate the critical strength of the repulsive interaction 
. Two 0~ states are found to have their critical strengths very close to 1, 
they are considered as candidates of new narrow resonances or loosely bound 
states . 

1, INTRODUCTION 

Since positrons can be easily created in the processes of high energy 
collisions, the molecules formed by electrons and positrons are believed to 
exist in nature. Half a century ago the existence of the positronium molecule, 
namely the dipositronium Ps2 . has already been predicted 1 and has been 
firstly calculated via a variational procedure 2 . Since the lifetime of the ground 
state of Ps2 is very short (it is only 0.906ns 3 ) due to the e~-e + annihilation, it 
has not yet been observed experimentally. In recent years this problem has 
attracted increasing attention following the ability to create cold positron 
beams tunable over a wide energy range 4 , and the increasing ability to carry 
out accurate theoretical calculations 3 ' 5 ~ 10 . In addition to the ground state, 
an angular momentum L=l and parity n = —1 bound excited state has also 
been predicted 10 . The other excited states are believed to be resonances. 

Of course any resonance would collapse via at least one open chan- 
nels. However, whether a state is allowed to enter into an open channel is not 
only determined by energy but also by symmetry. In fact, as we shall see, 
the states with a specific set of quantum numbers are allowed by symmetry 
to get access to only a few specific open channels. Thus, for a given state, it 
is important to make sure which open channel is the lowest one allowed by 
symmetry, this is called a symmetry-adapted lowest open channel (SLOC). 
If a state has its energy lower than the SLOC, it would be definitely bound 
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because it has no channel to collapse. If a state has its energy a little higher 
than the SLOC and if there is a barrier to hinder the wave function from 
leaking into the channel, it would appear as a narrow resonance. The barrier 
may have different origin, what is concerned here is the centrifugal barrier as 
discussed later. If the barrier does not exist or if the energy is much higher 
than the SLOC, the state would be a broad resonance. When a resonance has 
a very broad width, it is very difficult to be detected, and thus the existence 
of this resonance is meaningless. In other words, only the resonances with 
narrower widths are interesting. 

The main aim of this paper is the identification of the SLOC, this 
is the base for the further study of the dipositronium. Additionally, based 
on the SLOC and on other existing data, the existence of a few more bound 
states and resonances has been confirmed. Furthermore, a few candidates of 
bound states or resonances has been suggested based on an evaluation. The 
emphasis is placed on the qualitative aspect. 

2, SYMMETRY- ADAPTED LOWEST OPEN CHANNELS 

Let the particles 1 and 2 be electrons and 3 and 4 be positrons, let 
r a = r 2 -ri, r b = r 4 -r 3 , r c =|(r 4 +r 3 -r 2 -ri) , and 7^= | r ^- — r» | , The internal 
Hamiltonian of the dipositronium is (in a.u.) 

tf = -(v2 + vg + iv») + (£ + £- (i) 

which is invariant with respect to the O3 group and to the point group 
D 2d . 3 The latter is a subgroup of S 4 containing eight elements, namely 1, 
P12, P34, P12P34, P13P24, PUP23, p c (1324), and p c (1423) (where denotes 
an interchange and p c {ijkl) denotes a cyclic permutation). Hence the eigen- 
states can be classified according to the irreducible representations of the 3 
and D 2d , and thereby can be denoted as L n (fi), where \i denotes a repre- 
sentation of the T>2d group, pi — A±, A 2 , B±, B 2 or E. Let the wave function 
of a L n (fi) state be denoted as ^lu.^- From the knowledge of the D 2d group, 
we have 2 

^Ln M = *Ln> (if pi ^ E), or = - ^ Ln ^ (if fi = E) (2) 

and 

P13P24 ^Lii/i = ^Ln M (if H = Ax or B x ), or = - ^ LU ^ (if pi = A 2 or B 2 ) 

(3) 

Let the spins of the two electrons 1 and 2 be coupled to Si, those of 3 
and 4 be coupled to S 2 . Obviously, eq.(2) implies that the states with pi ^ E 
have Si — S 2 . 

On the other hand, let denotes an eigenstate of the positronium 
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Ps ( the ground state has (nl) = (10) ). For the Ps 2 ^Ps + Ps dissociation 
channel (where the Ps may be excited), The channel wave function \l/ c ^ can 
be written as 

V ch = (l+(-l) 5l Pi 2 )(l+(-l) 52 P34)[(^(ri3)^(r 2 4))zJ^(r 13 , 2 4)]L (4) 
where fi R is the wave function of relative motion and l R is the relative 

angular momentum. I and V are coupled to l Q , together with l R they are 

coupled to L. From (4) we have 

Pl3P24^= (-l) Z+ ^l + (-l) 52 Pl 2 )(l + (-l) 5l P34)[($nKri3)$nT(r 24 ))/J/ fl (ri3, 2 4)]L 

(5) 

Furthermore, in the special case of (nl) = (n'V) , eq.(4) can be rewritten 

as 

^ ch = (1 + (-l) Sl+52 ^ + ^){[($ n Kri3)^(r24)) i J/ fl (r 1 3, 24 )]L 

+(-l) 5l [(^(r23)^(ri4))/ /^(r 23 ,i4)]L} (6) 
Evidently, if a L n (ji) state is allowed to enter into a specific Ps + Ps 
channel, the associated ^/ c h must have the same (L, IT, //) symmetry. 

As examples, let us first study the 0~(Ai) and 0~(B 1 ) states. From 
eq.(3) we have 

Pl3P24^ch = ^ch (7.1) 

Since the \i = A\ or B\ states have S\ = S' 2 , eq.(5) can be rewritten as 

Pl3P2^ch = (-lY +l '^ch (7.2) 

Comparing eq.(7.1) with (7.2), we have 
(-1) ,+, ' = 1 (8.1) 

Furthermore, since we have meanwhile Si = S' 2 , the factor (l+(— l) Sl+S2+l ° +l R) 
in eq.(6) is nonzero only if l R + l is even. Thus we have 
(-1)Io+Ir =1? if (nl) = [n'V) (8.2) 

Besides, it can be proved that none of the partial waves of a 0~ state can 
be zero. The proof is as follows. Let the angular momenta corresponding to 
the three Jacobian vectors r , r&, and r c be denoted as l±, / 2 , and I3 . For CT 
states, any relative partial wave, say l±, can not be zero. If l\ were zero, 
I2 must be equal to l 3 to assure L=0. However, if l 2 — I3 , the parity must 
be even. Therefore l\ (or any relative angular momentum) is not allowed to 
be zero in 0~ states. Evidently, the proof is valid for any set of Jacobian 
coordinates. Thus we have 

l^L 0, V ^ 0, and l R ^ (8.3) 

Definitely, a Ps+Ps channel is accessible to the CT(Ax) or CT(£>i) states 
only if (8.1) to (8.3) are satisfied. Due to these constraints, the lowest Ps 
+ Ps channel is the one having (nl) = (n'V) = (21) , and l R = 1 . This 
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channel is labeled as (21) 2 (21), where p denotes the relative p-wave of the 
two excited positronium, the associated threshold energy is -0.1250 (a.u. are 
used in this paper). 

For the 0~(A 2 ) and 0~(B 2 ) states, the constraint eq.(8.1) should 
be changed to (-l) l+l '= -1 ( thus l R is even), while the constraints (8.2) 
and (8.3) remain unchanged. Accordingly the lowest Ps+Ps channel is the 
(21)-(32) channel at -0.090278. From a similar deduction, we know that 
the lowest Ps+Ps channels for the 0~(E) states are the (21)-(32) channel 
and the (21)2(31) channel both at -0.090278. 

Let us inspect another 2-body channel, the Ps~ + e + (or Ps + + e~ ) 
channel. It is recalled that the Ps - has only one bound state with angular 
momentum zero. Thus, if a 0~(fi) state entered into this channel, Ir must be 
zero to assure L=0. However, Ir = is prohibited due to eq.(8.3). Therefore 
the Ps~ + e + channel is not accessible to the 0~(fi) states. 

Let us inspect the Ps +e + + e~ 3-body channel. Let / be the an- 
gular momentum of the Ps. Since / = is not allowed due to eq.(8.3), The 
positronium in the 3-body channel of 0~ states must be excited, therefore 
the lowest Ps +e + + e~ 3-body channel is at -0.0625 (in which the Ps has 
(nl) = (21)). 

The above analysis is straight forward to be generalized and thereby 
all the SLOC can be identified as listed in Table 1 ( the states having sufficient 
data from existing theoretical calculations are listed in this table) and Table 
2, ( otherwise). 

Table 1, The SLOC and the corresponding threshold energies (in a.u.) of 
the L u (fj,) states (L < 2) together with the eigenenergies from theoretical 
calculations ( a is from ref.3 , b is from ref.10, and c is from our evaluation) 
. The partial wave given in the third column is the lowest wave. 



L 11 




SLOC 


threshold energy 


eigenenergy 


0+ 


Ai 


(10)s(10) 


-0.50000 


-0.516003778", -0.509 c 


0+ 


Bi 


(10)s(10) 


-0.50000 


-0.4994428°, -0.485 c 


0+ 


A 2 


(10)2(21) 


-0.31250 


-0.3120805", -0.301 c 


0+ 


B 2 


(10)2(21) 


-0.31250 


-0.3144689°, -0.307 c 


0+ 


E 


(10)^(20), (10)2(21) 


-0.31250 


-0.3300469°, -0.322 c 


\~ 


B 2 


(10)^(21) 


-0.31250 


-0.334408", -0.326 c 



Table 2, The SLOC, the threshold energies, and the critical strength A G . 
The A that smaller than 0.9 is not listed. 
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L 11 




SLOC 


threshold energy 


A 


0~ 


Ai 


(21)2(21) 


-0.12500 


0.93 


0" 


Bi 


(21)2(21) 


-0.12500 


0.98 


0" 


A 2 


(21)4(32) 


-0.090278 


0.94 


0" 


B 2 


(21)4(32) 


-0.090278 




0" 


E 


(21)2(31), (21)4(32) 


-0.090278 


0.99 


L 11 




SLOC 


threshold energy 






1+ 


Ai 


(10)4(32) 


-0.27778 




1+ 


Bi 


(10)4(32) 


-0.27778 




1+ 


A 2 


(10)2(21) 


-0.31250 




1+ 


B 2 


(10)2(21) 


-0.31250 




1+ 


E 


(10)2(21) 


-0.31250 




\- 


Ai 


(10)2(20) 


-0.31250 


0.93 


\- 


Bi 


(10)2(20) 


-0.31250 




\- 


A 2 


(10)2(21) 


-0.31250 




\~ 


E 


(10)2(10) 


-0.50000 


0.90 


L 11 




SLOC 


threshold energy 


A 


2+ 


Ai 


(10)4(10) 


-0.50000 




2+ 


Bi 

J -'l 


(10)4(10) 


-0.50000 




2+ 


A 2 


(10)2(21) 


-0.31250 




2+ 


B 2 


(10)2(21) 


-0.31250 


0.94 


2+ 


E 


(10)4(20), (10)2(21) 


-0.31250 


0.93 


2" 


Ai 


(10)2(32) 


-0.27778 




2" 


Bi 


(10)2(32) 


-0.27778 




2" 


A 2 


(10)^(21) 


-0.31250 




2" 


B 2 


(10)^(21) 


-0.31250 




2~ 


E 


(10)4(21) 


-0.31250 





3, DISCUSSION OF THE SPECTRUM 



(A) Confirmation of new bound states and resonance 
Among the states of Table 1, the O^A^ and 1~(B 2 ) have already 
been pointed out that they are bound. 3 ' 10 For the + (B 2 ) and + (E) states, 
they are tentatively classified as resonant singlet and triplet spin state in the 
first paper of ref. 7, and are believed to be metastable or resonance states 
in ref. 3. Now, it is clear that these two states have their energies lower than 
their thresholds, thus they are definitely bound. 

For the + (A 2 ) state, its energy is only a little higher than the thresh- 
old. Furthermore, the lowest partial wave for the relative motion of the two 
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dissociating positroniums is p— wave as shown in Table 1. Thus a centrifugal 
barrier exists to hinder the wave function from leaking out. It is recalled that 
the height of the centrifugal barrier of the Ps + Ps channel is ln(lR+l)/(2r^ h ), 
where r c h is roughly the sum of the radii of the two projectiles ( a positro- 
nium in the ground state and the other one in the (21) excited state), and 
therefore would be in the order of 10 . Accordingly the height of the barrier 
would be in the order of 0.01 for p— wave. On the other hand, the energy of 
the + (A2) is only higher than the threshold by 0.0004. 3 Thus the barrier is 
sufficiently high to hinder the wave function, and therefore we believe that 
the width of the + (A 2 ) resonance is very narrow. For the + (Bi) state, its 
energy is also only a little higher than the threshold. However, the lowest 
partial wave for the dissociation is s— wave, thus centrifugal barrier does not 
exist, and therefore the + (Bi) is a broad resonance difficult to detect. 

(B) Candidates of new bound states and resonances 

In order to evaluate the eigenenergies of the states in Table 2, a 
set of basis functions is introduced for the diagonalization of the Hamil- 
tonian. It is well known that the harmonic oscillator (h.o.) states are in 
general not appropriate for Coulomb systems, the main shortcoming is their 
inappropriate asymptotic behavior. However, if several set of h.o. states 
with different widths are used together, this shortcoming can be remarkably 
cured 11 . For example, let ip^f be the eigenstate of a single-particle pure 
harmonic oscillation with the width ujk , eigenenergy (2n + I + 3/2)hujK and 
angular momentum ITi . When five sets of ip^f are used together (i.e., ujk 
has five choices), and n < 3 is assumed for diagonalizing the Hamiltonian 
of a hydrogen atom, the resultant lowest energies for / = 0, l,and 2 states 
are -0.499987, -0.124998, and -0.055554, respectively, to be compared with 
the exact values -0.5, -0.125, and -0.055556. Therefore, we believed that, if 
only bound states are taken into account, and if only the qualitative aspect 
is concerned, several sets of h.o. states together can be used for our limited 
purpose. 

In the follows, the following seven sets of basis functions 
K K = te(r a )« 2 (r b )^ 3 (r c )]J i (9) 

are used together, where u>k = &i to uj 7 , ujijuj^x = 2.5, and (_i) z i +i 2+'3 
= n . Based on , the basis functions of a specific representation jj, can 
be induced by the following idempotents of the D 2c z group 3 . 

e M = |(1+P13P24)(1+P1 2 )(1+P34) (10.1) 
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e^ = f(l-p 13 p24)(l-Pi 2 )(l-P34) (10.2) 

e Bi = |(l+p 13 p 2 4)(l-Pl 2 )(l -P34) (10.3) 
e^ 2 = 1(1-P13P24)(1+P12)(1+P34) (10.4) 

e Eu = l(l-p 12 p 34 +p 13 p 24 -p 14 p 23 ) (10.5) 
e E 2 2 = |(l-p 12 p 34 -P13P24+P14P23) (10.6) 

Besides, a subsidiary procedure is needed to extract a new set {^i} from 
the old sets, so that the are orthonormalized. The total number of ^ 
used for the diagonalization is determined by N max which is the maximum of 
the sum 2(ni+n2+n 3 )+/ 1 +/ 2 + ^3 • As examples, when N max =23, the number 
of for the 0~(Ai) and 0~(Bi) states are 6655 and 7824 respectively. 

Evidently, the above basis functions are designed for bound states 
but not for resonances. However, most states in Table 2 are resonances (as 
we shall see). In order to avoid such an awkward situation, we introduce 
an adjustable parameter A in the repulsive interactions as A(^ + ^) ,while 
the attractive interactions remain unchanged. Such an adjustment does not 
change the threshold energies of the Ps+Ps channels. It is emphasized that 
the introduction of A does not at all alter the symmetry of the Hamiltonian, 
it remains to be invariant with respect to O3 and ■ Therefore the nature 
of the problem is not altered. Of course, what we really concern is the case 
with A close to I. 

Let the lowest eigenenergy of a given (LU/i) symmetry be denoted 
as E(\). Then, our procedure is firstly to choose a A smaller than one so 
that E(\) is lower than the threshold of the SLOC. For an example, for the 
CT(£>i) state, when A = 0.95, we have E(X) = —0.1282 which is lower than 
the SLOC at -0.1250. Secondly, A is increased step by step (in each step A 
is increased by 0.01). When A is equal to a value A c so that E(\ ) is lower 
than while E(X + 0.01) is higher than the threshold, then the procedure 
stops and A G is call a critical strength. In this way mainly bound states are 
concerned in the calculation. Evidently, if a state has a A G much larger than 
one, it is deeply bound. If A D is a little larger than one, it is just bound and 
a little lower than the dissociation threshold. If A D is a little smaller than 
one, it is a resonance a little higher than the threshold. If A G is remarkably 
smaller than one, the state is a high-lying resonance with a broad width, and 
we will neglect it. 

To show the convergency of our calculation, the eigenenergies of se- 
lected states are listed at the last column of Table 1 (with A = 1). Besides, 
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two series of eigenenergies are given in Table 3 in accord with N max . Ow- 
ing to the difficulty in calculation, the basis functions with N max ^23 are not 
adopted. The above results exhibit that the speed of convergency is not 
good, thus the calculated energies would not be useful if we want to know 
the exact locations of the levels (this is an effort beyond the aim and scope 
of this paper). However, if we just want to know which levels are relatively 
closer to their SLOC, then our results are useful. 

Table 3, The calculated energies E(\) (in a.u.) of the lowest 0~(Ai) state 
with A = 1 and the lowest 0~(Bi) state with A = 0.95 when N max is given. 



N r 



0-(B 1 



15 



-0.1088 



-0.1258 



The numerica 



17 



-0.1110 



-0.1267 



results o 



19 



-0.1128 



-0.1274 



21 



-0.1142 



-0.1279 



23 



-0.1153 



-0.1283 



A c are summarized in Table 2. In this table 
most A are remarkably smaller than one, the associated states are broad 
resonances and will not be further discussed. However, there are two and 
only two states, namely the 0~(E) and 0~(B 1 ) , have their A very close to 
one. They are distinguished from the others. Since, as mentioned before, all 
the calculated energies are a little higher than the actual values, the actual 
critical strength X ac tuai will be a little larger than the A listed in the table. 
Thus, the \ ac tuai of the 0~(E) and 0~(Bi) will be in fact very close to one 
or even larger. Furthermore, the open channels of these two states have Ir 
to be p— wave or higher, thus a centrifugal barrier exists. Therefore, among 
all the other states, these two states are candidates, they are either narrow 
resonances or bound states. If they are resonances, the 0~ (E) would emerge 
when two excited positroniums collide with each other (via the (21) 2 (31) 
or (21) -(32) channel). Similarly, the 0~(Bi) would emerge in the collision 
via the (21) -(21) channel. Once they are created, they are free from direct 
annihilation. 12 They would collapse either via their channel of formation, or 
would be transformed to a lower broad resonance via an El transition (the 
0~(E) would then become a 1 + (E), and the 0~(£>i) become a 1 + (A 2 )), and 
then collapse via the (10) 2 (21) channel. Anyway, during the decay of these 
resonances, low-energy photons (from the above El transition and from the 
transition of an excited positronium to its ground state) can be detected. 
These low-energy photons would help the identification of these states. 
4, FINAL REMARKS 

In conclusion, the main outcome of this paper is the identification of 
the SLOC of dipositronium. The identification is proved to be very useful to 
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the analysis of the existing data. Specifically, in addition to the two bound 
states that have been found previously, two more + bound states have been 
identified, and a + narrow resonances has been predicted. Furthermore, two 

states are found to have their energies very close to their SLOC (while 
the others are remarkably higher). Besides, their SLOC contain centrifugal 
barrier. Thus they are either narrow resonances or loosely bound states. 
Since our calculation is not accurate enough, what they really are remains 
to be identified . Undoubtedly, the final identification of these two states is 
an attractive theoretical topic. In the experimental aspect, the search of the 
+ (A 2 ) resonance via the collision taking place in the (10) 2 (21) channel is 
firstly recommended, because its width is very narrow. 

As a summary, a primary spectrum of dipositronium is proposed as 
shown in Table 4, where at least four bound states have been identified. 

Table 4, Bound states (b) and narrow resonances (r) of dipositronium 
together with their energies (in a.u.). The state lying upper is higher in 
energy. A is a small quantity (positive or negative) in the order of 0.001. 
The values of the energies of the + and 1~ states come from the literatures 
cited in Table 1. 





energy 




0~(E) 


-0,090278+A 


r or b 


O-(Bi) 


-0.1250+A 


r or b 


+ (A 2 ) 


-0.3120805 


r 


+ (B 2 ) 


-0.3144689 


b 


0+(E) 


-0.3300469 


b 


l-(5 2 ) 


-0.334408 


b 


o+(ao 


-0.516003778 


b 



The above procedure is mainly based on symmetry consideration, the 
way of analysis can be generalized to study other molecules with symmetries 
other than O3 and D 2 d- In particular, the identification of the SLOC is 
important to the study of resonances of various molecules. 
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